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Abstract By resorting to some results in quantum field theories with spontaneous break- 
down of symmetry we show that an explanation based on microscopic dynamics can be given 
of the fact that topological defect formation is observed during the process of non-equilibrium 
phase transitions characterized by a non-zero order parameter. We show that the Nambu- 
Goldstone particle acquires an effective non-zero mass due to boundary (finite volume) effects 
and this is related with the size of the defect. We also relate such volume effects with tem- 
perature effects. 



PACS: 03.70.+k, 11.30.Qc, 95.50.Eb 



Much attention is currently devoted to the problem of topological defect formation during 
the process of symmetry breaking phase transitions JT) . In such a process it may happen that 
a region, surrounded by ordered domains, remains trapped in the "normal" or symmetric 
state. This occurrence manifests as a topological defect. Topological defects are observed as 
macroscopic extended objects with classical behavior, e.g. vortices in superconductors and 
superfluids, magnetic domain walls in ferromagnets, dislocations, grain boundaries, point 
defects in crystals. In cosmology, topological defects, such as cosmic strings, may have been 
playing a role in the phase transition processes in the early Universe ||. The Kibble-Zurek 
scenario H Q provides the phenomenological understanding of the defect formation in phase 
transitions. There is a surprising analogy between defect formation in solid state physics and 
in high energy physics and cosmology ||]. For an interesting table of analogies see ||. As 
an example, we just mention the analogy between vortex in superfluids and global cosmic 
strings. Studying the physics of defect formation in condensed matter physics may be then 
helpful in the understanding of possible scenarios in the early Universe cosmology jfl]. The 
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analysis of the formation of defects in phase transitions thus becomes a "diagnostic tool" M 
in the study of non-equilibrium symmetry breaking processes in a wide range of energy scale. 

In this paper, by limiting ourselves to the class of non-equilibrium symmetry breaking 
phase transitions characterized by an order parameter, we discuss some of the questions 
which arise in the study of defect formation. For example, which one is the connection 
between defect formation and symmetry breaking, i.e., more precisely, why are extended 
objects with topological singularity observed only in systems showing some sort of ordered 
patterns? Why defect formation is observed during the processes of phase transitions? Why 
are the features of the defect formation so general, i.e. why are they shared by quite different 
systems, from condensed matter to cosmology? How to explain in a unified theoretical 
scheme the macroscopic behavior of extended objects and their interaction with quanta? 
How temperature and finite volume (the system boundaries) effects fit in such a theoretical 
frame? 

Of course, there is a large body of literature in quantum field theory (QFT) on the 
dynamics of defect formation and some of the above questions have been, at least partially, 
answered. By resorting indeed to some of such works, in the present paper we reach the 
following main results: 

i) we show that a theoretical explanation based on microscopic dynamics can be given of 
the fact that topological defect formation is observed during the process of non-equilibrium 
symmetry breaking phase transitions (limited to the case of the existence of the order pa- 
rameter and in the presence of a gauge field). 

ii) we provide an explicit proof that the Nambu-Goldstone (NG) particle acquires an 
effective non-zero mass due to boundary (finite volume) effects, which is related with the size 
of the defect. We also relate such volume effects with temperature effects. 

Our derivation is based on symmetry properties of the dynamics at the basic level of the 
quantum fields. As far as we know, in the literature the conclusions we reach are commonly 
accepted mostly on the basis of phenomenological arguments. Also, our conclusions are 
model independent, which may explain why some of the features of defect formation, within 
the limitations imposed to our analysis, are shared by many systems in a wide range of energy 
scale, independently of specific aspects of the system dynamics. Moreover, our results apply 
to the case of the Abelian U(l) symmetry (for vortices) as well as to the non-Abelian case of 
SO(3) and SU(2) symmetry (as for the monopole and the sphaleron case, respectively). As 
already mentioned, throughout this paper our discussion is limited to the case in which an 
order parameter exists and a gauge field is present. However, we also present some comments 
on boundary and temperature effects on the kink defect. 

We arrive at the result i) as a direct consequence of two known properties of gauge theories 
with spontaneous breakdown of symmetry (SBS). Both these properties have to do with the 
topological characterization of the non-homogeneous boson condensation in SBS theories. 
Although these properties are known since several years, nevertheless they have not been 
considered till now in relation with the problem of the appearance of topological defects in 
phase transitions. For a preliminary discussion see ref. [0]. For the reader convenience we 
report in the Appendices some details of the mathematical formalism to which we have to 
resort in our discussion. 

Let us start by recalling that in QFT the formation of extended objects is described by 
non-homogeneous boson condensation ||] - This condensation is formally obtained by 
translations of boson fields (not necessarily massless), say Xin{x) — » Xin(%) + f(x), with c- 
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number function f(x), satisfying the same field equation for Xin( x )- These translations are 
called boson transformations ||. The boson transformation function plays the role of a "form 
factor" and the extended object appears as the macroscopic envelop of the non- homogeneous 
boson condensate (localized over a finite domain). The topological charge of the extended 
object (the defect) arises from the topological singularity of the boson condensation function. 

We remark that, on the other hand, transitions between the system phases characterized 
by different ordered patterns in the ground state are induced by variations (gradients) of the 
NG boson condensate, i.e. of boson transformation function. Hence, in order to show that 
in such a process the conditions are met for the formation of topological defects, we need 

(a) to consider under which constraints the boson transformation function f(x) can carry 
a topological singularity and 

(b) to show that in the process of phase transitions under study these constraints are in 
fact satisfied. 

Let us first consider the point (a). We recall that in theories where boson condensa- 
tion can occur, topological defects are observed only in systems exhibiting ordered patterns 
(characterized by a non-zero order parameter). 

In fact one can show that boson transformation functions carrying topological singularities 
are allowed only for massless bosons ||, pX|] , such as NG bosons of SBS theories where ordered 
ground states appear. The proof of this goes as follows |l(J. 

Consider the boson transformation Xin( x ) ~^ Xin( x ) + f( x )- Let f(x) carry a topological 
singularity, which means that it is path-dependent: 



[d„,d„}f(x)^0 



for certain a , u , x . 



(1) 



We will see below that f is related with observables and therefore it is single- valued, 
i.e. [dp, d v \ d fl f(x) = 0. Recall that f{x) is solution of the Xin equation and suppose there is 
a non-zero mass term: (d 2 + m 2 )f(x) = 0. From the regularity of d fJi f{x) it follows that 
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d 2 + 



in- 



;d x Gl(x) 



(2) 



which leads to d 2 f(x) = 0, which in turn implies m = 0. Thus (fy is compatible only with 
massless Xin- This explains why topological defects are observed only in systems exhibiting 
ordered patterns, namely in the presence of NG boson condensation sustaining the long range 
ordering correlation. 

We remark that the above conclusion is not limited by dimensional considerations or by 
the Abelian or the non-Abelian nature of the theory symmetry group. It can be easily ex- 
tended to the non-relativistic regime (see [{Tcfl ) and it holds true in a theory with or without 
gauge fields. It applies to a full set of topologically non-trivial extended objects, ranging 
from topological line singularity, to surface singularity, point defects, grain boundaries and 
dislocation defects in crystals, SU(2)-triplet model and monopole singularity. The topologi- 
cal singularity and the topological charge of the related extended object can be completely 



characterized. For a detailed account see ref. [10], especially Chapter 10. 

Next, for the point (b), we consider the case in which a gauge field is present. We then 
observe that from Eq. (|2^) (Appendix B) we obtain the classical ground state current j M as 



Mx) = {0\j H ^x)\0) 



my 



a^{x) d ll f{x) 

eo 



(3) 
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where m v a fl (x) is the Meissner-like current and ^-d^fix) is the boson current. 

Note that the observable classical current is given in terms of gradients d u f of the boson 
condensation function. Eq. (||) shows that the macroscopic ground state effects, such as the 
classical field and the classical current, do not occur for regular f(x) (CrL, = 0). In fact, from 
( |28| ) we obtain d 2 a a (x) = for regular f(x), i.e. a a (x) = ■^d ll f(x), which in turn implies 
zero classical current (j„ = 0) and zero classical field {F„ v = d u a u — d u a^) (the Meissner-like 
current and the boson current cancel each other). 

The conclusion is that the vacuum current appears only when f(x) has topological singu- 
larities, which, as we have seen above (point (a)), is only compatible with the condensation 
of massless bosons, i.e. when SBS occurs. Again this conclusion is a general one, it does not 
depend on the specific model. It holds for Abelian as well as for non-Abelian groups, in the 
relativistic and in the non-relativistic case. 

Summing up, in a gauge theory, the symmetry breaking phase transitions characterized 
by macroscopic ground state effects, such as the vacuum current and field, (e.g. in super- 
conductors) can occur only when there are non-zero gradients of topologically non trivial 
condensation of NG bosons. But this means that the conditions for the formation of topo- 
logical defects are met and thus we see why topological defects are observed in the process of 
symmetry breaking phase transitions. We have thus obtained our result i). 

Notice that the appearance of space-time dependent order parameter v is no guarantee 
that persistent ground state currents (and fields) will exist: if / is a regular function, the 
space-time dependence of v can be gauged away by an appropriate gauge transformation. 
Which means that non-homogeneous condensation, and thus extended objects, cannot be 
obtained in such case. Also note that in a theory which has only global gauge invariance 
non-singular boson transformations of the NG fields can produce non-trivial physical effects 
(like linear flow in superfluidity). 

In conclusion, in a gauge theory with non-zero order parameter, in order to have topo- 
logical defects one needs boson condensation function with topological singularities and this 
is only allowed for NG boson condensation. These conditions are all met in the process of 
phase transition in such a theory. 

A further consequence of the above analysis is that, since the boson transformation with 
regular / does not affect observable quantities, the S matrix must be actually given by (cf. 
Appendix A) 

S = : S[p in , U? n - —d( X in ~ h n )] : • (4) 
my 

This is in fact independent of the boson transformation with regular /: 

S _> S ' =: S[p in , U? n - —d( X in - b in ) + ZP(a" - -d»f)] : (5) 

my eo 

since a n {x) = ■^d fl f(x) for regular /. However, S' ^ S for singular /: in such a case (||) 
shows that S' includes the interaction of the quanta U^ n and pi n with the classical field and 
current. This shows how it may happens that quanta interact and have effects on classically 
behaving macroscopic extended objects. 

For the point ii), let us consider now the effects of finite volume on homogeneous boson 
condensation. 

For large but finite volume we expect that the symmetry breakdown order parameter is 
constant "inside the bulk" far from the boundaries. However, "near" the boundaries, one 
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might expect "distortions" in the order parameter: "near" the system boundaries we may 
have non-homogeneous order parameter, v = v(x) (or even v — * 0). Such non-homogeneities 
in the boson condensation will "smooth out" in the V — > oo limit. 

From the Ward-Takahashi identities one obtains the following pole structure for the two- 
point function of the x( x ) ^eld fl3fl: 



ix( x )x(y)) = nm { 7 — / d p—^ — 5 h (continuum contributions) . (6) 

e-»o I (27r) 4 J p z — m x + iea x 

Z x and a x are renormalization constants. The space integration of (x( x )x(y)) then picks up 
the pole contribution at p 2 = 0, and leads to |l4], |l5| 

Z 



J x 



v 44> m x = , or v = m x 7^ , (7) 



where v denotes a convenient c-number [13]. Eq. (j7|) expresses the well known Goldstone 
theorem: if the symmetry is spontaneously broken {v 7^ 0), the NG massless mode exists, 
whose interpolating Heisenberg field is Xh{x)- Since it is massless it manifests as a long range 
correlation mode and it is thus responsible for the vacuum ordering, as already observed 
above. 

Restrict now the space integration of Eq. (|6|) over the finite (but large) volume V = n~ 3 
and use for each space component of p: 

5Jp) = — [" dxe ipx = — sin^ , (8) 
nV ' 2vr J-l -np n W 

v 

which, as well known, approaches 5(p) as n — > 0: lim v ^o5 v (p) = 5(p). Consider that 

lim^o J dp5 v (p) f(p) = /(0) = lim v ^ J dp5(p-n) f(p) . (9) 

Then, using S r] (p) ~ 5(p — n) for small n, one obtains 

v(y, e, n) = ieve'^ A x (e, n, Po = 0) , (10) 



with 



A x (e,r/,p = 0) 



" ' + (continuum contributions) 



u 2 + tea. 
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and io 2 p=ri = rf + m x . Thus, lim t ^olim v ^ov(y , e, n) 7^ only if m x = 0, otherwise v = 0. 
Note that the Goldstone theorem is thus recovered in the infinite volume limit (77 — > 0). 

On the other hand, if m x = and n is given a non-zero value (i.e. by reducing to a finite 
volume, i.e. in the presence of boundaries), then u> p=71 7^ and it acts as an "effective mass" 
for the x bosons. Then, in order to have the order parameter v different from zero e must be 
kept non-zero. 

In conclusion, near the boundaries (77 7^ 0) the NG bosons acquire an effective mass 
m e ff = u> p=r) . They will then propagate over a range of the order of £ = -, which is the 
linear size of the condensation domain, or, in the presence of topological singularity, the size 
of the topologically non-trivial condensation, namely of the extended object (the defect). 
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We stress that if n ^ then e must be non-zero in order to have the order parameter 
different from zero: v 7^ (at least locally). In such a case the symmetry breakdown is 
maintained because e / 0: e acts as the coupling with an external field (the pump) providing 
energy. Energy supply is required in order to condensate modes of non-zero lowest energy 
cjp =7? . In summary, boundary effects are in competition with the breakdown of symmetry 
||. They may preclude its occurrence or, if symmetry is already broken, they may reduce to 
zero the order parameter. 

Temperature may have similar effects on the order parameter (at Tc symmetry may be 
restored). Since the order parameter goes to zero when NG modes acquire non-zero effec- 
tive mass (unless, as seen, external energy is supplied), we may then represent the effect of 

thermalization in terms of finite volume effects and put, e.g., n oc \] ^ T ~t^ c ^ - In this way 
temperature fluctuations around Tq may produce fluctuations in the size £ of the condensed 
domain (the size of ordered domain and/or of the domain where non- homogeneous conden- 
sation occurs, namely of the defect). 

For example, suppose to start with T > Tc, but near to Tc- In the presence of an external 
driving field (e ^ 0), even before transition to fully ordered phase is achieved (as T — » Tq), 

one may have the formation of ordered domains of size £ oc {\f^^~')~ 1 ■ As far as n ^ 0, the 
ordered domains (and the topological defects) are unstable, they disappear as the external 
field coupling e — > 0. 

Of course, assuming T is lowered below Tq, if ordered domains are still present at T < Tc, 
they also disappear as e — > 0. The possibility to maintain such ordered domains below Tc 
depends on the speed at which T is lowered, compared to the speed at which the system is 
able to get homogeneously ordered. Notice that the speed at which T — * Tc is related to the 
speed at which 77—5-0. 

In the case of the kink solution, the boson transformation function fp(x) = const. e _W) (^) Xl 
plays the role of "form factor" for the soliton solution. The number of condensed bosons is 
proportional to \fp(x)\ 2 = e~ 2 ^ ^ (xi-a) ^ w j 1 j c j 1 [ s maximal near the kink center x\ = a and 
decreases over a size £g = j^hn ■ The mass po = (2\)zv(j3), with []ll]] v 2 {f3) = v 2 — 3(: p 2 : 
}o < v 2 , of the "constituent" fields p m (x) fixes the kink size £3 oc = , which thus 

increases as T increases: T — > Tc (say T ^ Tc but near Tc- At Tc, v((3c) = 0.). 

In the T — > limit the kink size is £n oc -7=- < = fn. For T different from zero, 

the thermal Bose condensate (: p 2 :}o develops, which acts as a potential term for the kink 
field. It is such a potential term which actually controls the "size" (and the number) of the 
kinks. Only far from the kink core, namely in the limit v(x, f3) — > const, the p in {x) field may 
be considered as a free field. The translation of the boson field by fp breaks the homogeneity 
of the otherwise constant in space order parameter v{(5). 

Let us close by observing that in the case of topologically non-trivial condensation at 
finite temperature the order parameter v(x,f3) provides a mapping between the variation 
domains of (x, (3) and the space of the unitarily inequivalent representations of the canonical 
commutation relations, i.e. the set of Hilbert spaces where the operator fields are realized 
for different values of the order parameter. This expresses the known fact that we have non- 
trivial homotopy mappings between the (x, (3) variability domain and the group manifold. As 
well known, in the vortex case one has the mapping -k of S , surrounding the r = singularity, 
to the group manifold of U(l) which is topologically characterized by the winding number 



6 



n € Z € tti(S 1 ). It is such a singularity which is carried by the boson condensation function 
of the NG modes. In the monopole case the mapping tt is the one of the sphere S 2 , 
surrounding the singularity r = 0, to SO(3)/SO(2) group manifold, with homotopy classes of 
^2(S 2 ) = Z. Same situation occurs in the sphaleron case [11], provided one replaces SO(3) 
and SO(2) with SU(2) and U(l), respectively. 

In conclusion, transitions between phases characterized by an order parameter imply 
"moving" over unitarily inequivalent representations, and this is achieved by gradients in NG 
boson condensation function. In the presence of a gauge field, macroscopic ground state field 
and currents can only be obtained by non-homogeneous NG boson condensation with topo- 
logical singularities. In turn, the occurrence of such topologically non-trivial condensation 
allows the formation of topological defects: this explains why topological defect formation is 
observed in such a kind of symmetry breaking phase transition processes ( "where the defects 
come from"). 

We have also seen that finite volume (the presence of boundaries) is related to the size 
of the ordered domains and of the defects and we have briefly discussed the correlation with 
temperature effects. This may justify the expectation relating the number of observed defects 
with the non-equilibrium dynamics of phase transitions. 

In the case of the kink there are no NG modes, nevertheless the topologically non-trivial 
kink solution requires the boson condensation function to carry divergence singularity (at 
spatial infinity). 

This work has been partially supported by INFN, INFM, MURST and the ESF Network 
on Topological defect formation in phase transitions. 
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Appendix A. The boson condensation and the dynamical rearrangement of symmetry 



In the standard Lehmann-Symanzik-Zimmermann (LSZ) formalism [16| of QFT, the dy- 
namics is given in terms of the interacting fields, say ipn(x), also called the Heisenberg fields. 
On the other hand, observables are described by asymptotic in- (or out-) operator fields 
(physical fields), say (fi n (x), which satisfy free field equations. The mechanisms described 
in this appendix are fully general, however, for simplicity we consider the Lagrangian den- 
sity C[4>h{x), (pH*(x), Ajj^ix)] for a complex scalar field 4>h{x) interacting with a gauge field 
Ahii(x)- We do not need to specify the detailed structure of the Lagrangian. We only require 
that such a Lagrangian be invariant under global and local U(l) gauge transformations (as, 
e.g., in the Higgs-Kibble model |17|, p|]): 



(f) H (x) 



4> h {x) -> e^^Hix) 



A Hfl (x) -» A Hfl (x) , (12) 

A Hfl (x) -» A HfjL {x) + d fJi X(x) , (13) 

respectively, with A(x) — > for \xq\ — > oo and/or |x| — > oo. We use the Lorentz gauge 
d^An^x) = and put ((>h{x) = [ipu{x) + ixn(x)]. We also assume that SBS can occur: 
(0\4>h(x)\0) = v 7^ 0, with v constant and put Ph{x) = iPh(x) — v. It is known [13 1 that the 
theory contains a massless negative norm field (ghost) bi n (x), the Nambu-Goldstone massless 
mode Xin( x ), and a massive vector field U^ n . The relation between the dynamics and the 
observable properties of the physical states is provided by the LSZ mapping (the dynamical 
map, or Haag expansion) between interacting fields and physical fields. In the present case 



one finds [13] the following LSZ maps: 
<j) H (x) =: exp 



i—Xinix) 
v 



?p Pin(x) + F[p in , Uf^, d(xin ~ h n )] 



(14) 



A»(x) = ZlU? n (x) + ^2Ldn in (x)+ : F»[p in ,U? n ,d( Xt 
eov 



(15) 



Z x , Z p and Z3 are the wave function renormalization constants. As usual, the colon symbol 
denotes normal ordering. The functionals F and F^ have to be determined within a specific 
choice for the Lagrangian. The S'-matrix is given by 5 = : S[pi ni U^ n ,d{xi n — bin)] '■■ It 
is important to notice that the dynamical mappings (0) and (ID (and the one for the S'- 
matrix) are weak equalities, i.e. they are equalities among matrix elements computed in the 
Fock space for the physical states. The free field equations are 



9 2 Xin(x) 







b in (x) = 
{d 2 + m v 2 )U? n (x) = 



(d 2 + m 2 )p in (x) 



d»U? n (x) 







(16) 
(17) 



with my 2 = ^-(e v) 2 . 

One can show that the local gauge transformations and the global phase transformations 
(cf. (0) and (|i~2|)) of the Heisenberg fields are induced by the in- field transformations (see 
@ and dlj)): 



Xin(x) -> Xin(x) + ^T-A(x) , 
Zx 



bin(x) -> b in (x) + ^rA( 



(18) 
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Pin(x) - Pin (x) , Ut n {x) - U? n (x) . (19) 

and by 

Xin(x) -> Xin(^) + , (20) 

6in(«) -> &inO*0 -» Pin(x) ,Ui n {x) -> Uj* n (x), (21) 

respectively. The square integrable function /(x) is required to be solution of d 2 f(x) = 



0. Note that ( p0| ) with /(x) = 1 (translation by a constant c-number) is not unitarily 
implementable. f(x) is introduced in order to make the generator of such a transformation 
well defined. The limit f(x) — > 1 (i.e. the infinite volume limit) is to be performed at the end 
of the computation. Notice that the in-field equations and the S matrix are invariant under 
the above in-field transformations (in the limit / — ► 1 ). 

The transformation ( p0|) with f(x) = 1 describes the homogeneous boson condensation. 
Once NG bosons are condensed in the ground state, long range correlation is established which 
manifests as the ordered pattern in the system ground state. This means that transitions 
between the system phases characterized by different ordered patterns in the ground state 
may be induced by the process of boson condensation. In other words, phase transitions 
(transitions among unitarily inequivalent Fock spaces) are induced by variations (gradients) 
of the NG boson condensate. 

The fact that the local and the global U(l) gauge transformations are induced by the 



group G' of transformations (18-21) is named the dynamical rearrangement of symmetry 
|], 10 1 : The dynamical rearrangement of symmetry expresses the consistency between the 
invariance of the Lagrangian under the symmetry transformations and the SBS condition 
(Q\4> H {x)\Q)=v^Q. 

Translations of boson fields (eqs. ( |l8| ) and (p0[)) are thus obtained as a consequence of 
SBS. As a general result one finds that G' is the group contraction of the symmetry group of 
the dynamics (U(l) in our present case; more generally, if the symmetry group is SU(n) or 
SO(n), the group contraction G' is EU(n-l) or E(n), respectively) p9| . Notice that G' is the 
transformation group relevant to the phase transitions process. 



Appendix B. Macroscopic field and currents generated by the boson condensation 

In the framework of the model in Appendix A, the Maxwell equations are 

-d 2 A» H {x) = f H {x) - d^B(x) , (22) 

with 

B(x) = ^Ws) - Xin(x)] , d 2 B{x) = , (23) 

and with jfj(x) = 5£(x)/5Ah^(x). We may require that the current jjj is the only source of 
the gauge field A^j in any observable process. This is obtained by imposing the condition: 
p (b\d^B(x)\a) p = 0, i.e. 

(-d 2 ) p (b\A^ H (x)\a) p = P (b\f H (x)\a) p , (24) 
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where we use A°^(x) = A^(x) — e$v : d^bi n (x) :. \a) p and \b) p denote two generic physical 
states. The condition p {b\d tl B{x)\a)p = leads to the Gupta-Bleuler-like condition 



o . 



(25) 



with \\n an< ^ the positive-frequency parts of the corresponding fields. Xin and bi n do 
not participate to any observable reaction. However, we stress that the NG bosons do not 
disappear from the theory: their condensation in the vacuum can have observable effects. 
Notice that Eq.(24) are the classical Maxwell equations. 

We now remark that the boson transformation must be also compatible with the physical 
state condition (|25|). Under the boson transformation Xin(x) — * Xin(x) + B changes 



as 



7^ 



B(x) 



B{x) 



e v 2 



/(*) • 



(26) 



Eq. (24) is then violated. In order to restore it, the shift in B must be compensated by 
means of the transformation on Ui n : 



Ut n (x)^U? n (x) + Z 3 ->a»(x) 



o , 



(27) 



with a convenient c-number function a^(x). The dynamical maps of the various Heisenberg 
operators are not affected by (|27]) provided 



(d 2 +my)aJx) 



m 



V 
eo 



daf{x) • 



(28) 



This is the Maxwell equation for the vector potential a^ [[0|, 12]. A similar result may be 
obtained in more complex models with non-Abelian symmetry groups [10, |i"l| ]. 
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